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Abstract 

The dimensionally regularized massless double box Feynman diagram 
with powers of propagators equal to one, one leg off the mass shell, 
i.e. with non-zero q 2 = pf, and three legs on shell, p 2 = 0, i = 2,3,4, 
is analytically calculated for general values of q 2 and the Mandelstam 
variables s and t. An explicit result is expressed through (generalized) 
polylogarithms, up to the fourth order, dependent on rational combi- 
nations of q 2 , s and t, and a one-dimensional integral with a simple 
integrand consisting of logarithms and dilogarithms. 
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1 Introduction 



Massless four point Feynman diagrams contribute to many important physical ampli- 
tudes. They are much more complicated than two- and three-point diagrams because 
depend on many parameters: the Mandelstam variables s and t and the values of the 
external momenta squared, p 2 , i = 1,2,3,4. In the most general case, when all the 
legs are off the mass shell, p 2 7^ 0, there exists an explicit analytical result [|TJ for the 
master (i.e. with powers of the propagators equal to one) double box diagram (see 
Fig. 1) strictly in four dimensions. Still no similar results are available for pure off 




shell four point diagrams with ultraviolet, infrared and/or collinear divergences. 

In the opposite case, when all the end-points are on shell, i.e. for p 2 = 0, i — 
1,2,3,4, the problem of the analytical evaluation of such diagrams, in expansion in 
e = (4 — d)/2 in the framework of dimensional regularization with the space-time 
dimension d cis du rc gularization parameter, was completely solved during last year in 
H, U |J. Among intermediate situations, when some legs are on shell and the rest 
of them off shell, the case of one leg off shell, q 2 = p 2 ^ and three legs on shell is 
very important because of the relevance to the process e + e~ — > 3jets (see, e.g., ||). 
The purpose of this paper is to analytically evaluate the master double box diagram 
of such type, as a function of q 2 ,s and t, and thereby demonstrate that the NNLO 
analytical calculations for this process are indeed possible. 

One of the ways to evaluate the four point diagrams with one leg off shell is 
to expand them in the limit q 2 —>■ and compute as many terms of the resulting 
expansion as possible. We explain how to do this, following the strategy of regions 
J?]. H , in the next section and present the leading power term in this expansion which 
provides a very non-trivial check of the subsequent analytical result. 

To analytically evaluate the considered diagram we straightforwardly apply the 
method of ref. |§: we start from the alpha-representation of the double box and, after 
expanding some of the involved functions in Mellin-Barnes (MB) integrals, arrive at 
a six-fold MB integral representation with gamma functions in the integrand. Then 
we use a standard procedure of taking residues and shifting contours to resolve the 
structure of singularities in the parameter of dimensional regularization, e. This 
procedure leads to the appearance of multiple terms where Laurent expansion in e 
becomes possible. Resulting integrals in all the MB parameters but the last two are 
evaluated explicitly in gamma functions and their derivatives. The last two-fold MB 
integral is evaluated by closing an initial integration contour in the complex plane to 
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the right, with an explicit summation of the corresponding series. A final result is 
expressed through (generalized) polylogarithms dependent on rational combinations 
of q 2 , s and t and a one-dimensional integral with a simple integrand consisting of 
logarithms and dilogarithms. 

2 Expansion in the limit q 2 —> 

The dimensionally regularized master massless double box Feynman integral with one 
leg off shell, q 2 = p\ ^ 0, and three legs on shell, p 2 = 0, i = 2, 3, 4, can be written as 



F{s,t,q 2 ;e) = j J 



d d kd d l 



(k 2 + 2 Pl k + q 2 )(k 2 - 2p 2 k)k 2 (k - I) 2 
1 



x ■ 



(l 2 + 2p 1 l + q 2 ){l 2 -2p 2 l){l+p 1 +p- i ) 2 ' 

where s = {p\ +P2) 2 , t = {p\ +P3) 2 , and k and I are respectively loop momenta of 
the left and the right box. Usual prescriptions, k 2 = k 2 + zO, s = s + iQ, etc are 
implied. To expand the given diagram in the limit q 2 — > one can apply the so-called 
strategy of regions J7|, § based on the analysis of various regions in the space of the 
loop integration momenta, Taylor expanding the integrand in the parameters that are 
considered small in the given region and extending resulting integrations to the whole 
integration domain in the loop momenta. When applying this strategy all integrals 
without scale are by definition put to zero. 

Let us choose, for convenience, the external momenta as follows: 

q 2 

Pl=Pl ~Q2 P2 > Vi=V2, Pi,2 = (=FQ/2,0,0,<5/2), 

where s = —Q 2 . The given limit \q 2 \ \s\,\t\ is closely related to the Sudakov limit 
so that it is reasonable to consider each loop momentum to be one of the following 
types: 



hard (h): k ~ Q ~ y—t 



1-collinear (lc): k + ~ q 2 /Q, k- ~ Q , k ~ y—q 2 



2-collinear (2c): k + ~ Q, k^ ~ q 2 IQ , k ~ y—q 2 ■ 

Here k± = k ± k 3 , k— (k%, k^- We mean by k ~ Q, etc. that any component of k^ 
is of order Q. 

It turns out that the (h-h), (lc-h) and (lc-lc) are the only non-zero contributions 
to the leading power behaviour in the limit q 2 — > 0. Any term originating from the 
(h-h) contribution is given by the expansion of the integrand in Taylor series in q 2 and 
expressed through on-shell double boxes in shifted dimensions and can be analytically 
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evaluated by the algorithm presented in Q]. The (lc-lc) contribution is obtained by 
expanding propagators number 2, 4 and 7 in a special way. In particular, propagators 
number 2 and 4 are expanded, respectively, in I 2 and k 2 . (See || for instructive 2-loop 
examples of expansions in limits of the Sudakov type.) 

The (lc-h) and (lc-lc) contributions are evaluated with the help of a two-fold 
(respectively, one-fold) MB representation. Still this program of the evaluation of a 
large number of terms of the expansion looks very complicated because one needs, 
for phenomenological reasons, the values of q 2 greater than s and t so that a reliable 
summation of a resulting series, using Pade approximants, requires the knowledge of 
at least first 20-30 terms. Such a great number of terms can be hardly evaluated 
since a lot of irreducible structures appear. This asymptotic expansion is however 
very useful for comparison with the explicit result derived below. 

The leading power terms of the asymptotic expansion calculated in expansion in 
e, up to a finite part, are 

. 2 



F{s,t,q 2 ;e) 



\2+2e(_ 



t) 



E 

i=0 



+ 0(q 2 \n 3 (q 2 /s)) + 0{e) 



(2) 



where X = q 2 / s , Y = t/s and 



9i{X,Y) 
9s{X,Y) 

92(X,Y) 
9i{X,Y) 



-2(lnX-lnF) 

117T 2 



12 



+ 31nX lnY- -ln z Y 
2 



21nFLi 2 (-F) -2Li 3 (-Y) 



— ln 3 F + ln 2 F ln(l + Y) + vr 2 
6 



-ln 3 X--ln 2 XlnF--lnX In 2 Y 

3 2 2 

19 



lnX lnF + ln(l + r) 

6 



+ 



49C(3) 
6 



9o{X,Y) = 
+21nFLi 3 



26 Li 4 ( 

Y 
1 + Y 



■Y) -2S 2 , 2 (-Y) -2(lnX + 61ny + ln(l + F))Li 3 
+ (In 2 Y + 2 In X In Y + An 2 ) Li 2 (-Y) 



•Y) 



--ln 4 X + - ln 3 XlnF + - ln 2 Xln 2 F- -lnXln 3 F + - ln 4 F 

2 2 4 2 8 



+ ln(l + F) 



lnXln 2 F- -ln 3 F + - In 2 F ln(l + Y) 

3 2 



i 



InF ln 2 (l + Y) 



2 7 2^ 

— In 2 X - - lnX In Y + — In 2 Y + lnX ln(l + Y) - 2 In Y ln(l + Y) 

3 3 6 



+-ln 2 (l + F) 



+ C(3) 

Here Li a (z) is the polylogarithm 

S a ,b(z) 



19 34 

— lnX lnF + 21n(l + F) 

3 3 



+ 



83tt 4 
180" 



and 



1)!6! Jo 



dt 



(3) 



(4) 



the generalized polylogarithm 



3 From alpha parameters through MB represen- 
tation to analytical result 

The alpha representation of the double box looks like: 

2 P oo p oo 

F(s,t,g 2 ;e) = -r(3 + 2e) (m d/2 ) / de*i . . . / da 7 5 (^2 a i ~ l) D 1+3e A~ 3 ~ 2e , 

(5) 

where 

D = (ai + a 2 + a 7 )(a 3 + a 4 + a 5 ) + a e (ai + a 2 + a 3 + a 4 + a 5 + a 7 ) , (6) 
A = [aia 2 (a 3 + a 4 + a 5 ) + a 3 a 4 (ai + a 2 + a 7 ) + a 6 (a;i + a 3 )(a 2 + c*4)](— s) 

+a 5 a e a 7 (-t) + a 5 [(ai + a 3 )a 6 + a 3 (ai + a 2 + a 7 )](-g 2 ) . (7) 

As it is well-known, one can choose a sum of an arbitrary subset of , % = 1, . . . , 7 
in the argument of the delta function in (||), and we use the same choice as in 0. 

Starting from (|5|) we perform the same change of variables as in 0] and apply 
seven times the MB representation 



/ dw———T(is + w)T(-w) 



(X + Y) v T(v) 2ni J- 

in order to separate terms in the functions involved to make possible an explicit 
parametric integration. The two extra MB integrations arise form the extra term 
with q 2 . After such integrations we are left with a 7- fold MB integral of a ratio of 
gamma functions. Fortunately, one of the integrations can be explicitly taken using 
the first Barnes lemma and we arrive at the following nice 6-fold MB integral: 

(™ d/2 Y 1 r 

F(s,t, q 2 ; e) = — — — — — — / dvdwdw 2 dw 3 dzdz 




r(-l -3e)(-s) 3 + 2e (2ni) e 
xT(l + u>)r(l + v + w)T(-v)T(-w)T(l - w 3 + v) 
xr(tw 2 )r(-l -2e-w- w 2 )T(w 3 - v)V(-l - 2e - w - w 3 ) 
T(l - w 2 + zi)r(l - w 3 + Zi)V(e + w + w 2 + w 3 - zx)T(-zi) 
T(l + w + w 2 + iw 3 )r(— 1 - 4e — w — w 2 — w 3 )T(l - w 3 ) 
xT(l - e + z)T(2 + 2e + w + w 2 + z~ z 1 )T(2 + 2e + w + w 3 + z-zi] 

T(— 2 — 3e — w — w 2 — w 3 + z\ — z)Y{z\ — z) 
X r(3 + 2e + w + z) ' 



(9) 



It differs from its analog for q 2 = by the additional integration in v. This variable 
enters only four gamma functions in the integrand. The integral is evaluated in 
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expansion in e, up to a finite part, by resolving singularities in e absolutely by the 
same strategy as in the case q 2 = . Note that the infrared and collinear poles are 
a little bit softer than in the pure on-shell case, the integration variable v playing the 
role of an infrared regulator. The two key gamma functions that are responsible for 
the generation of poles in e are the same as in the previous case: 

r(e + w + w 2 + u>3 - zi) r(— 2 - 3e — w — w 2 — w 3 + z\ — z) . 

The labeling of resulting terms is therefore the same: the initial integral is decomposed 
as J = Jqo + Joi + Jio + Ju, etc (Only arguments of some gamma functions are shifted 
by v.) The applied strategy makes it possible to perform all the integrations apart 
from the last two, in v and w. We obtain four groups of terms with 26 terms in each 
group: the terms without MB integration, with MB integration in v or w and, finally, 
with a two- fold integration in v and w. The one-fold integrals are explicitly evaluated 



by closing contour and summing up series, using formulae from ] I 



The contribution of the resulting two-fold MB integral takes the form 



2 



2[m d l 2 ) 1 rdvdwfq 2 \ u ft 



-s 3 (2%i) 2 



2 



I 5 I 

T(l + v + w)r(-v-w) (- -7E-21n(-s) 



,e l + w 1 + v + w 

+ip(l + v) - 2ip(-v - w) - 3ip(-w) + 2^(1 + w) + -0(1 + v + w)) 

-T(l + v)T(-v)T(l + w)T(-w)} . (10) 

The integration contours are straight lines along imaginary axes with —1 <Ret>, 
Reu>,Ret> + w < 0. By closing contours it is possible to convert this integral into a 
two-fold series where each term is identified as a derivative of the Appell function F 2 
in parameters, up to the third order. The 1/e part is then explicitly summed up with 
a result in terms of polylogarithms. (In fact, it is proportional to the e part of the 
master one- loop box.) 

The so obtained result can be transformed into a one- dimensional integral with a 
simple integrand. To present the final result let us turn to the variables x = s/q 2 and 
y = t/q 2 keeping in mind typical phenomeno logical values of the involved parameters 
relevant to the process e + e~ — > 3jets: 

F (s, t, q\ e) = I _ . * £ + 0(e) . (11) 



-sH{-q 2 



i=0 



We obtain 



f 4 {x,y) = -1, (12) 
f 3 {x,y) = 2(\nx + \ny), (13) 
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f 2 {x,y) = 3Li 2 (x)+ Li 2 (y)-2(lnx + lny) 5 



+3 ln(l — x) lnx + ln(l — y) lny — 



5tt^ 
12 ' 



(14) 



= 2 

— lnx Li 2 



Li 3 

y 



1-x-yJ \l-x-y j 



Li 3 



-xy 



1 — x — y 



1 — x 



In |/ Li 2 



x 



21n(l-x-y) 



x 



- (ln 2 (l — x — y) + 7r 2 ) + ln(l — x) lnx + ln(l — y) lny — lnx lny 

+3 Li 3 (x) - 8 Li 3 (y) + 4 Li 3 (^) - 2 Li 3 ^^-) 

4 2 

-(3mx + 41ny)Li 2 (x) + 31nyLi 2 (y) + -In 3 a; - -ln 3 (l - x) + ln 2 (l - x) lnx 

o o 

9 7T 2 4 1 

--ln(l - x) ln 2 x + — (5 lnx -41n(l - x)) + -ln 3 y + -ln 3 (l - y) 
A o o 



7T 



-21n 2 (l -y) lny-ln(l - y) In 2 y + — (5 lny + 2 ln(l - yj) 



6 



+4 lnx lny (lnx — ln(l — x) + lny) + 



25C(3) 
6 



(15) 



The e° part involves a one-dimensional integral: 

fo(x,y) 

jo 
An 

[ln(l - yz) (ln(l - z) ln(l -yz)-2 Li 2 (z)) 



y dzjz 1 ln(l - z)(41n 2 (l - x - yz) - ln 2 (l - y - xz)) 
Ay 



1 — x — yz 

-2(ln(l - 2) - In z) Li 2 (-(1 - x - yz)/x)\ 

— [ln(l - xz) (3 ln 2 (l -z)-6 ln(l - z) ln(l - xz) + 2 Li 2 (2 

+2(6 ln(l - z) - In z) Li 2 (-(1 - y - xz)/y)}} 

x \ „ „. . ( —x \ „ T . ( xy \ 



-5 Li 4 (x) + 14 Li 4 



-2Li 4 



1 — x 



6 Li 4 



+8Li 4 — + 24Li 4 (y) - 2Li 4 (1 - y) + 8Li 4 

vi-x-yy V 1 -?/ 



(l-x)(l-y)y 
- 2Li 4 



y 



1 — X 
1 — x — y 
1 — x 



-8Li 4 (l -x) -8Li 4 ( ) -20Li 4 (- — - — - ) + 10 Li 4 

\l-x-yj v l ~v J 

-3S 2y2 (x) - 8S 2 , 2 (y) - QS 2 , 2 (j^j + (2 In y - 2 In x - 3 ln(l - x)) Li 3 (x) 
+2(16 ln(l -y) - 11 lny - 2 ln(l - x - y) + lnx) Li 3 (j—^j 



(81ny + 21n(l - x) + 31nx) Li 3 



-x 



1 — X 



xy 



-2(41n(l-y)-41ny-ln(l-x)+lnx)Li 3 , 

V(l -x)(l -y)J 

+ (14 ln(l - y) - 18 lny + 4 ln(l - a; - y))Li 3 ( ^"^ J + 2 lny Li 3 
+ (7 In y - 8 ln(l - y))Li 3 (y) + (8 ln(l - y) + In y + 2 In x)Li 3 
-4(21ny + 71n(l - or) + 21n(l - x - y) -81na;)Li 3 
-2(lny + 51n(l - x) + 81n(l - x - y) -71na;)Li 3 

1 



-xy 



I- x-y 



X 

-y 



(Li 2 (x)Y- Li 2 



^(Li 2 (y)) 2 + 4(Li 2u ^" X2 



1 — a; — y 
1/ 



+ 



ln 2 (l — a;) — 4 In 2 y + 2 lny (4 ln(l — x) — In x) — 2 ln(l — y) lna; — 3 ln(l — a;) In a; 
5tt 2 ~ 



7- 2 



Li 2 (x) + 12 ln 2 (l - y) + 15 In 2 y + 2 ln(l - y)(ln(l -x-y)+\nx 



+ 2 hl+ 3 
-91ny) + 21ny (ln(l - a; - y) - 41n(l - rr) + 4 In a) 



-2(ln 2 (l -x-y) + \n 2 x) Li 2 



x 



+ 



-41n 2 (l -y) - 11 In 2 y 



+ 

-81n(l -x)(ln(l - x-y) - 2 lna;) + 2 ln(l - y) lnx + ln^x- 
+ [ln 2 (l -a) -41n 2 (l - y) - 8 In 2 y + 2 lny (4 ln(l - a;) - 
+21n(l - y){A\ny - In a;) - 21n(l - a;) In a; + 21n 2 xl Li 2 





L12 (y) 


y) — 8 ln(l — a; — y) lna; 


2tt 2 1 

a; 

3 




3 In a;) 




xy 


-y)) 


(l-x)(l 





+2 ln 4 (l - x - y) + - ln 3 (l - x - y) [2 lny - 3 ln(l - y) - 9 lna; - 11 ln(l - a;)] 

+ ln 2 (l - a; - y) [tt 2 - 3 ln 2 (l - y) + 3 In 2 y + 6 ln 2 (l - a;) - ln(l - y)(lny - 10 lna;) 

+4 ln(l — x) lna; — 2 In 2 x — lny (5 ln(l — x) + lna;) H — ln(l — a; — y) 71n 3 (l — y) 

3 

-2 ln(l - y)(57T 2 + 12 In 2 y) + 7tt 2 ln(l - x) + tt 2 lna; + 6 In 2 y lna; - 4 ln 3 (l - a;) 
+ 151n 2 (l -y)(lny-21nx) -211n 2 (l - a;) lna; + 31n(l - x) In 2 a; 
+ In y (2tt 2 + 9 In 2 x + 15 In 2 (1 - x) - 6 ln(l - x) In x)" 
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5 2 23 17 7 

— ln 4 (l — x) In 4 x H ln 3 (l — x) lnx ln 2 (l — x) In 2 x H — ln(l — x) In 3 x 
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7T 2 2 

-— (31n 2 (l -x) - 101n(l - a;) In a; + 5 In 2 a;) - ln 4 (l - y) - -ln 4 y 

19 2 

- — ln 3 (l - y) lny + 5 ln 2 (l - y) In 2 y + - ln(l - y) In 3 y 

6 3 



7T 

~6~ 



9 ln 2 (l — y) — ln(l — y) lny — 5 In 2 y 



H — ln(l — x) ln(l — y)(ln 2 (l — x) — 4 ln 2 (l — y)) (In 2 x + In 2 y) lnx lny 

3 3 

+31n 3 (l -y)lnx + ln 2 (l - y) [lny(41n(l - x) + lnx) -21n(l - x)h\x] 

+— lny — ln 3 (l — x) — 9 ln 2 (l — x) lnx — 12 lny In 2 x 

o 

+6 ln(l — x) lnx (2 lny + 3 lnx)] — ln(l — y) 8 In 2 y ln(l — x) 



+ ln(l — x) lnx (ln(l 



x 



2 lnx) + lny (—8 ln 2 (l — x) + 61n(l — x) lnx + In 2 x) 



71 



H [In y (4 ln(l — x) — 5 lnx) — ln(l — y) lnx] 

3 



+C(3) 



25 



12(ln(l - x - y) - ln(l - y)) + 131n(l - x) (lnx + lny) 

3 



+ 



237T 4 

180" 



(16) 



One may hope that the one-dimensional integral that is left can also be evaluated in 
terms of polylogarithms. To do this it is necessary to complete the table of integrals 
derived in ||. 

This result is in agreement with the leading power behaviour when q 2 
When performing this comparison it is reasonable to start with (|TUD, take minus 
residue at v = (the first pole of T(— v)), integrate in w by closing the contour to the 
right, and take into account the three other contributions (without MB integration, 



and with integration in v or w) that were not presented above. Eqs. (|T2|--[16[) also 
agree with results based on numerical integration in the space of alpha parameters 
Ul2] (where the 1% accuracy for the 1/e and e° parts is guaranteed). 
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